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explain that (i) at low energy this theory does not exactly recover general relativity 
but instead mimic general relativity plus dark matter; that (ii) higher spatial curvature 
terms allow bouncing and cyclic universes as regular solutions; and that (iii) the 
anisotropic scaling with the dynamical critical exponent z = 3 solves the horizon 
O ■ problem and leads to scale-invariant cosmological perturbations even without inflation. 

We also comment on issues related to an extra scalar degree of freedom called scalar 
graviton. In particular, for spherically-symmetric, static, vacuum configurations we 
prove non-perturbative continuity of the A — > 1 + limit, where A is a parameter in the 
\jj kinetic action and general relativity has the value A = 1. We also derive the condition 

^v ■ under which linear instability of the scalar graviton does not show up. 
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1. Introduction 

One of the biggest difficulties in attempts toward the theory of quantum gravity is the 
fact that general relativity is non-renormalizable. This would imply loss of theoretical 
control and predictability at high energies. In January 2009, Hofava proposed a 
new theory of gravity to evade this difficulty by invoking a Lifshitz-type anisotropic 
scaling at high energy pp. This theory, often called Hofava-Lifshitz gravity, is power- 
counting renormalizable and is expected to be renormalizable and unitary. Having a 
new candidate theory for quantum gravity, it is important to investigate its cosmological 
implications. 

There are a number of interesting cosmological implications of Hofava-Lifshitz 
gravity. For example, higher spatial curvature terms lead to regular bounce solutions in 
the early universe [21 [3] . Higher curvature terms might also make the flatness problem 
milder [3]. The anisotropic scaling with z = 3 solves the horizon problem and leads to 
scale- invariant cosmological perturbations without inflation [5]. The anisotropic scaling 
provides a new mechanism for generation of primordial magnetic seed field [6], and 
also modifies the spectrum of gravitational wave background via a peculiar scaling of 
radiation energy density [Zj. In parity- violating version of the theory, circularly polarized 
gravitational waves can also be generated in the early universe jH]- The lack of local 
Hamiltonian constraint leads to dark matter as an integration "constant" [91 [TO]. 



The purpose of this article is to review basic construction of the theory and some 
of its cosmological implications. In Sec. [2] we explain basics of Hofava-Lifshitz gravity, 
such as power-counting argument, symmetry, basic quantities, action and equations of 
motion. In Sec. |3] we comment on issues related to an extra scalar degree of freedom 
called scalar graviton, and consider the limit in which general relativity is supposed to be 
recovered. We explicitly see the known result that the naive metric perturbation breaks 
down in this limit for the scalar graviton. However, this does not necessarily imply the 
loss of predictability. Indeed, for spherically-symmetric, static, vacuum configurations 
we shall show that the limit is non-perturbatively continuous. This result might 
correspond to what is called Vainshtein mechanism [TT] in theories of massive gravity [12] 
and suggest that the extra scalar degree of freedom might safely decouple from the rest 
of the world in the limit. In Sec. H]we shall review some of cosmological implications: the 
dark matter as an integration "constant" , bouncing and cyclic universes and generation 
of scale-invariant cosmological perturbation without inflation. Finally, Sec. [5] is devoted 
to a summary of this article and some discussions. 

2. Hofava-Lifshitz gravity 

2.1. Preliminaries 

2.1.1. Power- counting Let us begin with heuristically explaining the usual power- 
counting argument in field theory. As the simplest example, let us consider a scalar 
field with the canonical kinetic term: 

- fdtd 3 x<p 2 , (1) 

where an overdot represents a time derivative. The scaling dimension of the scalar field 
4> is determined by demanding that the kinetic term be invariant under the scaling 

t -^-bt, x -► bx, 4> -¥ b~ s (j), (2) 

where b is an arbitrary number and s is the scaling dimension to be determined. The 
invariance of the kinetic term under the scaling leads to the condition 

1 + 3 - 2 - 2s = 0, (3) 

where 1 comes from dt, 3 from d 3 x, —2 from two time derivatives and —2s from two 
0's. Thus, we obtain s — 1. In other words, the scalar field scales like energy. With this 
scaling in mind, it is easy to see that an n-th order interaction term behaves as 

' dtd z x<p n oc E- (1+3 - ns) , (4) 



where E is the energy scale of the system of interest. Here, the minus sign in the 
exponent comes from — 1 in E — > b~ x E, 1 in the parentheses comes from dt, 3 from 
d 3 x and — ns from <fi n . Now, it is expected that we have a good theoretical control of 
ultraviolet (UV), i.e. high E, behaviors if the exponent is non-positive. Since s = 1, 
this condition leads to n < 4. This is the power-counting renormalizability condition. 



We are interested in gravity. Unfortunately, Einstein gravity is not power-counting 
renormalizable. This is because the curvature is a highly nonlinear functional of the 
metric and there are graviton interaction terms with n higher than 4. The non- 
renormalizability is one of difficulties in attempts to quantize general relativity. 

2.1.2. Abandoning Lorentz symmetry As already stated, Hofava-Lifshitz gravity is 
power-counting renormalizable. How does it evade the above argument? The basic idea 
is very simple but potentially dangerous: abandoning Lorentz symmetry and invoking 
a different kind of scaling in the UV. The scaling invoked here, often called anisotropic 
scaling or Lifshitz scaling, is 

t ->■ b z t, x -> bx, <f) ->■ b~ s (j), (5) 

where z is a number called dynamical critical exponent. 

Let us now see how the power-counting argument changes if the scaling is 
anisotropic as in (J5]). Invariance of the canonical kinetic term ([T]) under this scaling 
leads to 

z + 3 - 1z - Is = 0, (6) 

where z comes from dt, 3 from d 3 x, —2z from two time derivatives and —2s from two 
0's. Then we obtain 

* = ¥• < 7 » 

This of course recovers the previous result s — 1 for z — 1. What is interesting here is 
that s = if z = 3. This implies that, if z = 3, the amplitude of quantum fluctuations of 
4> does not change as the energy scale of the system changes. The n-th order interaction 
term behaves as 

dtd 3 x<f) n oc E- {z+3 - ns)/z , (8) 

where —1/z in the exponent comes from — z in E — > b~ z E, z in the parentheses comes 
from dt, 3 from d 3 x and — ns from <fi n . For z = 3 (and thus s — 0), the exponent 
is negative for any n and, therefore, any nonlinear interactions are power-counting 
renormalizable. For z > 3, the theory is power-counting super-renormalizable. 

From the above consideration, it is expected that gravity may become 
renormalizable if the anisotropic scaling with z > 3 is realized in the UV. 

2.1.3. Scalar field action We would like to realize the anisotropic scaling with z > 3 
in the UV to construct renormalizable nonlinear theories. On the other hand, in order 
to recover the Lorentz invariance in the infrared (IR), we would like to realize the usual 
scaling with z = 1 at low energy. A simple example is a scalar field with the following 
free-part action: 

I free = \j dtd 3 X (</> 2 + <f)0<f)) , (9) 



where 

A 3 kA 2 o a o 

° = W-Tp +c * h -< < 10) 

M is the energy scale corresponding to the transition from the z = 1 scaling to the 
z = 3 scaling, k is a constant, c^ is the sound speed, i.e. the limit of speed in the IR, 
m^ is the mass of the field, and A is the Laplacian in the 3-dimensional space [j. 

In the UV, the sixth-order spatial derivative term dominates over lower-order terms 
and balances with the time kinetic term which includes two time derivatives. This 
naturally leads to the z = 3 scaling. On the other hand, in the IR, the second-order 
spatial derivative term and the mass term are dominant and, thus, the z = 1 scaling is 
realized. In this way, it is possible to realize the z = 3 scaling in the UV and the z — 1 
scaling in the IR. 

However, one must be aware that all "constants" in the action are subject to running 
under the renormalization group (RG) flow. Of course, the sound speed is not an 
exception. If we consider many fields then the sound speed for each field should run 
under the RG flow [15) . We need a mechanism or symmetry to make sound speeds of 
different species to be essentially the same at low energies. More generally speaking, we 
need a mechanism or symmetry to suppress Lorentz violating operators at low energies. 
Perhaps, embedding the theory into a larger theory is necessary. One such possibility 
is related to supersymmetry [16] . 

2.2. Symmetry 

As explained in the previous subsection, the way the power-counting renormalizability 
is achieved is to violate the Lorentz invariance and to invoke the anisotropic scaling with 
the dynamical critical exponent z > 3. Since the Lorentz invariance is not respected, 
we treat the time coordinate t and the spatial coordinates x 1 (i = 1,2,3) separately. 

The fundamental symmetry of the theory is the invariance under space-independent 
time reparametrization and time-dependent spatial diffeomorphism: 

t^t\t), x->tf(t,x). (11) 

The time-dependent spatial diffeomorphism allows an arbitrary change of spatial 
coordinates on each constant time surface. However, the time reparametrization here is 
not allowed to depend on spatial coordinates. As a result, unlike general relativity, in 
Hof ava-Lifshitz gravity the foliation of spacetime by constant time hypersurfaces is not 
just a choice of coordinates but is a physical entity. Indeed, the foliation is preserved 
by the symmetry transformation (ITT1) . For this reason, the map ([lip is called foliation 
preserving diffeomorphism. 

In addition to the foliation preserving diffeomorphism invariance, we assume that 
the theory is invariant under the spatial parity x — > —x [17] and the time reflection 
t-> -t. 

| If photons have this kind of dispersion relation with k = 0(1), c<f, = 1 and m<p = then experiments 
such as Fermi GBM/LAT [13] and MAGIC [14, set a lower bound on M as M > 10 u GeV. 



Finally, in order to render the theory power-counting renormalizable, we would like 
to realize the anisotropic scaling with z > 3 at high energy. In the present article, for 
concreteness, we mainly focus on the case with z = 3. 

2. 3. Basic quantities and projectability condition 
Basic quantities of Hof ava-Lifshitz gravity are 

lapse : N(t), shift : N\t,x), 3d metric : ^ (£,£), (12) 

from which we can construct a four- dimensional spacetime metric of the ADM form as 

ds 2 = -N 2 dt 2 + g^dx 1 + N i dt){dx j + N j dt). (13) 

While the shift N l and the 3d metric gij depend on both the time coordinate t and the 
spatial coordinates x, the lapse N is assumed to be a function of the time only. This 
condition on the lapse is called the projectability condition. 

The projectability condition stems from the foliation preserving diffeomorphism. 
The lapse represents a gauge freedom associated with the space-independent time 
reparametrization t — > t'(t) and, thus, it is fairly natural to restrict it to be space- 
independent l§|. Of course, the projectability condition is compatible with the foliation 
preserving diffeomorphism. The transformation of the basic quantities ( 1121) under the 
infinitesimal foliation preserving diffeomorphism, 

St = f{t), 5x i = Cit ) x) ) (14) 

is defined as follows. 

6N = d t (Nf), 

5N l =d t (N i f) + d t ? + £tN i , 

S(N { ) = dt(Nif) + gydt? + CsN i7 

$9ij = fdt9ij + C^g ij} (15) 

where Ni = gijNK Note that 5N is independent of spatial coordinates since / and N 
are functions of time only. Thus the projectability condition is compatible with the 
foliation preserving diffeomorphism: the foliation preserving diffeomorphism maps a 
space-independent iV to a space-independent N. 

The equation of motion for the lapse corresponds to the generator of the time 
reparametrization and is called the Hamiltonian constraint. Since the lapse is 
independent of spatial coordinates, its variations are also space- independent. This 
means that the Hamiltonian constraint in Hofava-Lifshitz gravity is not a local equation 
but an equation integrated over a whole space. In subsection 14.11 we shall discuss 
cosmo logical implication of the global nature of the Hamiltonian constraint. 

§ Abandoning the projectability condition leads to phenomenological obstacles [18] and theoretical 
inconsistency [19]. On the other hand, the criticisms made in |181fT9] do not apply if the projectability 
condition is respected. 



Kij = ^ (k - D t N 3 - D 3 N t ) . (17) 



2.4. Action 

The theory should respect the foliation preserving diffeomorphism. We can then use the 
following ingredients in the action: 

Ndt, y/gd^x, gij, A, R ih (16) 

where g is the determinant of g^, Di is the 3- dimensional covariant derivative compatible 
with Qij and Rij is the Ricci tensor of g^. Note that the Ricci tensor includes all 
information about the Riemann tensor since Weyl tensor identically vanishes in 3- 
dimensions. 

2.4-1 ■ The UV action We should include time-derivative of the 3-dimensional metric 
in the action in order to make the metric dynamical. However, g io is not covariant under 
the spatial diffeomorphism and, therefore, g^ should appear in the action as a part of 
the covariant quantity called extrinsic curvature, 

1 
2iV 

The extrinsic curvature transforms as a second-rank symmetric tensor under the spatial 
diffeomorphism and as a scalar under the time reparametrization. The time kinetic term 
for the metric is obtained by squaring the extrinsic curvature and properly contracting 
indices. There are two ways to contract indices: 

hvn = ^Q J Ndt^d 3 x {K^K l3 - \K 2 ) , (18) 

where G and A are constants, and K = K\. In general relativity, A is fixed to 1 
because of higher symmetry. On the other hand, in Hofava-Lifshitz gravity, any value 
of A is compatible with the foliation preserving diffeomorphism invariance and thus A 
is not fixed. We shall not include terms including derivatives of the extrinsic curvature 
in the action. This is consistent if the theory without those higher derivative terms 
is renormalizable: those terms would be non-renormalizable and thus would not be 
generated by quantum correction. For the same reason, we shall not include terms 
cubic or higher order in the extrinsic curvature. 

Invariant terms made of time derivatives of the shift would inevitably include second 
or higher time derivatives of the spatial metric. For the reason explained above, we shall 
not include those higher derivative terms in the action. Time derivative of the lapse 
corresponds to the connection in 1-dimension spanned by the time but the curvature in 
1-dimension is always zero. Thus, there is no invariant term made of time derivatives 
of the lapse. Of course, the spatial derivative of the lapse vanishes because of the 
projectability condition. 

Since terms in the kinetic action (1181) include two time derivatives, we should include 
terms with six spatial derivatives in order to realize the z = 3 scaling in the UV. (For a 
general choice of z (> 3) in the UV, we should include terms with 2z spatial derivatives.) 



The foliation preserving diffeomorphism invariance allows five such terms in the action, 
I z=3 = f Ndt^d 3 x [dDiRjkD 1 ^ 

+c 2 D i RD i R + csRitfRi + c±RR\R) + c 5 R 3 ] , (19) 

where q (i — 1, • • • , 5) are constants. Note that the other possible term DiRjkD^ R kl is a 
linear combination of the above terms up to total derivative and, thus, does not have to 
be included explicitly. We do not include terms with more than six spatial derivatives 
since they would be non-renormalizable and thus would not be generated by quantum 
corrections if the theory without them is renormalizable. 

2.4-2. Relevant deformations in the IR In the IR, terms with less number of spatial 
derivatives in the action become important. There are two independent terms with four 
spatial derivatives 



I z=2 = / Ndt^/gd 3 x [cqRIR] + c 7 R 2 \ , (20) 

one term with two spatial derivatives 

I z=1 = c 8 JNdt^d 3 xR, (21) 

and a constant 

h=o = c 9 J Ndt^d 3 x , (22) 

where q (i — 6, • • • , 9) are constants. 

We have written down all possible terms consistent with the symmetry of the theory 
except for terms involving more than two time derivatives and terms with more than 
six spatial derivatives. As already stated, those higher- derivative terms excluded in the 
above construction would be non-renormalizable and, thus, would not be generated by 
quantum corrections if the theory without them is renormalizable. The theory defined 
in this way is power-counting renormalizable and, thus, expected to be renormalizable 
although renormalizability beyond the power-counting argument has not been proved. 
Also, the theory is expected to be unitary since the action does not include more than 
two time derivatives. Note that the constants G, A and q (i — 1, • • • , 9) are subject to 
running under the RG flow. 

2-4-3. IR action with z = 1 In the UV the theory naturally exhibits the z = 3 scaling 
as the second time derivative terms Ikin and the sixth spatial derivative terms I z= s 
balance with each other. 

On the other hand, in the IR the forth and sixth spatial derivative terms, I z= 2 
and I z= 3, are unimportant. We therefore have the following action describing the IR 
behavior of the theory: 



Ml, 



IlR — Ikin + Iz=l + h=0 

j Ndt^/gd 3 x (K ij Kij -XK 2 + R- 2A) , (23) 
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where Mj, z = 1/(8ttG) and A = — 87rGcg, and we have set 167rGcg to unity by rescaling 
of the time coordinate. This IR action naturally exhibits the z = 1 scaling. Moreover, 
the action looks identical to the Einstein-Hilbert action in the ADM formalism if A — > 1. 
There are however two important differences: (i) A does not have to be 1 and is subject 
to running under the RG flow; (ii) the projectability condition restricts the lapse A" to 
be a function of the time only. Regarding (i), the RG flow of the theory has not been 
investigated and, thus, we do not know whether A = 1 is an IR fixed point of the RG 
flow or not. On the other hand, we shall discuss cosmological implication of the point 
(ii) in subsection 14.11 

2. 5. Equations of motion 

Adding the matter action I m , the total action is 

/ =I g + I m , (24) 

M 2 r 
I g = — ^ / Ndt^gd z x{K l3 K %j -XK 2 + A + R + L z>1 ), (25) 

M 2 

— ^L z>1 = (dDiR^D'R^ + ^DiRD'R + c 3 ^^^ 

I J 

+ cJiR\R) + c 5 R 3 ) + (ceRiR) + c 7 R 2 ). (26) 

Here, we have rescaled the time coordinate to set 167rGcs to unity. Note that not only 
the gravitational action I g but also the matter action I m should be invariant under the 
foliation-preserving diffeomorphism. 

By variation of the action with respect to the lapse N(t), we obtain the Hamiltonian 
constraint 



H 9± + H m± = 0, (27) 



where 



and 



W -Jd xH g± , H mX = - — 



H g± = --£ = / d 3 xU g± , H m± = --£, (2f 



M 2 
U g± = ^^(KV PlJ -A-R-L z>1 ), PtJ =K t3 -XKg ir (29) 

Variation with respect to the shift N l (t, x) leads to the momentum constraint 

T-Lgi + T-Lmi — 0, (30) 

where 

U gi = -^- = -M 2 Pl JgDi Pij , n mi = - 5 J^. (31) 

Note that the gravitational part of the momentum constraint is determined solely by the 
kinetic terms and thus is totally insensitive to the structure of higher spatial curvature 
terms. In particular, for A = 1 the momentum constraint agrees with that in general 
relativity. 



For comparison, let us consider the case in which the matter sector recovers 
spacetime diffeomorphism invariance. In this case it makes sense to define the stress- 
energy tensor T^ v of matter and then 



H ml _ = I d 3 xy/g T fMV n tJ "n v , H mi = —T^ri 



V9 



where 



1 



n^dx* = -Ndt, n^ = -Jd t - N l di). 



(32) 



(33) 



As in general relativity, the gravitational action can be written as the sum of kinetic 
terms and constraints up to boundary terms: 



dtd 3 x 



where 7r Jjf is momentum conjugate to g^ given by 



7T 



i.l 



6l n 



n ij d t gij - N'Hgi - / dtNH gL + (boundary terms), (34) 



(35) 



M 2 Pl ^gp l - 



p'- 



g ik g jl Pki- 



5(dt9ij) 

The Hamiltonian corresponding to the time t is the sum of constraints and boundary 
terms as 

H g [dt] = NH g± + J d^xWUgi + (boundary terms). (36) 

Finally, by variation with respect to gij(t,x), we obtain the dynamical equation 

gij ' ^"rnij U, yo I ) 

where 



£ ■■ 



2 5I„ 



9ik9jV 



£, 



2 SI„ 



N^g Sg kl ' -"" ~ 9ik9jl N^g Sgl ~ ^ ' (38) 

Note that the matter sector (as well as the gravity sector) should be invariant under 
spatial diffeomorphism (as a part of the foliation preserving diffeomorphism) and thus 
it makes sense to define T t j in general. The explicit expression for £ gi j is given by 



£ •• 



M 2 Pl 



~(dt - N k D k ) Pij + ^( Plk D 3 N k +p 3k D t N k ) 



■Kpij + 2Kfp kj + - gij K kl p kl + -A gij - dj 



£ 



z>liji 



(39) 



where £ z> uj is the contribution from L z>1 and Gij is Einstein tensor of gfy. 

The invariance of I a under the infinitesimal transformation ( TT5|) leads to the 
following conservation equations, where a represents g or m. 



= Nd t H a± + / d 3 x 
= hd t - NW, 



N l d t U ai + -Ny/gS^dtQij 



N 



2hA +K n 

V9J 



1 He 



V9 N ^g 



DtN 3 - D'Eaij. 



(40) 
(41) 
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3. Scalar graviton and the A — > 1 + limit 

3.1. Propagating degrees of freedom 

In order to identify propagating degrees of freedom, let us consider linear perturbations 
around a flat background without matter. We can decompose the perturbation into 
scalar, vector and tensor parts, according to the transformation under infinitesimal 
spatial diffeomorphism. Thus, we have 

N =1 + A, Ni = diB + m, 

ga = (1 + 2C)6ij + 2didjh L + dihj + djhi + hij, (42) 

where n, and hi are transverse and hy is transverse traceless: S^difij = 5 l: *dihj = 0, 
S^dihjk = and S^hij = 0. Also, A depends only on t because of the projectability 
condition. By fixing the gauge degrees of freedom /(£) and £;(£, x) as / = — / Adt and 
£i = —ditiL — hi, the gauge transformation ffT5|) leads to 

N = l, Ni = diB+7H, g ij = (l + 2£)8 ij + h ij . (43) 

In this gauge, the momentum constraint ( |30l) without matter is 



(3A - 1)C - (A - l)d 2 B] + l -d 2 n t = 0, (44) 



leading to 



B = T^ih n4 = ' (45) 

where d 2 = b % ^d{dj is the spatial Laplacian. We do not have to solve the Hamiltonian 
constraint since it is an equation integrated over a whole space and thus does not reduce 
the number of local physical degrees of freedom. The scalar physical degree of freedom 
( is often called scalar graviton while the tensor perturbation h^ represents the two 
physical degrees of freedom of usual tensor graviton. 

The time kinetic term expanded up to quadratic order is 

"3A-l- 9 1 



i-kin 



m 2 pi J dtd 3 x -j^C 2 + g FWhijhu 



+ 0(e J ), (46) 



where we have introduced a small expansion parameter e and considered ( and h^ as 
0(e). In order to avoid ghost instability, thus A must be either larger than 1 or smaller 
than 1/3 PQ. Since general relativity has A = 1 and we would like to recover something 
similar to general relativity in the IR, we should consider the regime A > 1. Although 
the RG flow of the theory has not yet been analyzed, a hope is that the RG flow may 
have a UV fixed point at A = +oo and an IR fixed point at A = 1 + so that A runs 
from +oo in the UV to 1 + in the IR. 
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3.2. Dispersion relation 



Expanding the potential terms up to the second order and adding them to ( 146 
obtain 

~3A-1 



we 



M; 



Pl 



dtd x 



A-l 



C + CO s C + 



h ik v l (hijhu + h tj o t h kl 



where 



a 



A-l 



'^1 

Mt 



K*A< 



A 



O t 



Mt 



k,A 2 
Mi 



A. 



3A- 1 \M* M 2 S 

Here, we have introduced M Pl = l/(8nG), set 167rGcg to unity by rescaling of the time 
coordinate, set A = in order to allow the flat spacetime as a consistent background, 
and defined M si and k s j as 

AC 4 = - 2(3ci + 8c 2 )Mpf , 

k s M; 2 = - 2(3c 6 + 8c 7 )M; 

Thus the dispersion relation is 
A-l / k e 



+0(e 3 ),(47) 



(4? 



'p; > 



Mr 4 

K< M- 



-2ciMpi 
-2c 6 Mp, 2 



p/ 



(49) 



CJ 



3A-1 VM S 4 
for scalar graviton, and 
k 6 



CK/atXt , O 



UJ 



K t k 



Mt 



M? 



(50) 



(51) 



for tensor graviton. 

As we have already seen, the absence of ghost requires A > 1. The dispersion 
relation (150]) then implies that the scalar graviton is unstable for k lower than ~ 
M J20j EU [22] and that the time scale of this linear instability is 



tL ~k\ 



3A-1 



(52) 



A-l 

As we shall see in subsection I4.1[ the lack of local Hamiltonian constraint leads to "dark 
matter as an integration constant", a non-dynamical component which behaves like 
pressure-less dust. As in the standard cold dark matter (CDM) scenario, the dust-like 
component exhibits Jeans instability and forms large-scale structures in the universe. 
The timescale of Jeans instability is 
M Pl 



t, 



Vp 



(53) 



where p is the energy density at the position of interest. Note that this instability is 
necessary for structure formation if we consider the dust-like component as an alternative 
to CDM. Thus, as far as 



t L > tj, 



(54) 
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the linear instability of the scalar graviton does not show up. Also, the linear instability 
is tamed by Hubble friction if 

t L > H~\ (55) 

where H is the Hubble expansion rate at the time of interest. If either (154]) or (155]) is 
satisfied then the linear instability of the scalar graviton does not show up [23]. For 
length scales shorter than ~ 0.01mm, we do not experimentally know how gravity 
behaves and, thus, the linear instability at shorter length scales would not contradict 
with any experiments. Also, modes with k higher than ~ M s are stable, provided that 
3ci + 8c 2 < 0. 

In summary, the condition under which linear instability of the scalar graviton does 
not show up is 



0< 



A-l 

3A-1 



< max 



IP 

k 2 ' 



l$l 



for H < k < min 



M e 



0.01mm 



(56) 



where we have introduced Newton potential $ by M Pl k 2 § ~ —p. Note that A is subject 
to running under the RG flow and thus should depend on k, H and $ in general. 
Therefore, the condition ( 156]) should be considered as a phenomenological constraint on 
properties of the RG flow. 



3.3. Breakdown of metric perturbation 

Basically, the condition ( 156]) says that A (> 1) must be sufficiently close to 1 at low 
energy, while A — 1 (> 0) can be of 0(1) or larger at high energy. In the following we 
shall show that a naive metric perturbation breaks down when A is close to 1. Non- 
perturbatively, however, the theory is described by a finite number of parameters, Mpi, 
A, A and q (i = 1, 2, • • • , 7) if renormalizable. 
A natural nonlinear extention of ( H3]) is 

,2C 



jV = l, Ni^diB + m, g tj 



where n^ is transverse and hij is transverse traceless: S^diUj 



(57) 

0, 5 lj dihj k = and 

S^hij = 0. We shall consider ( and h^ as O(e) and perform perturbative expansion 
with respect to e. 

In order to calculate the action up to cubic order, it suffices to solve the momentum 
constraint up to the first order. Thus, by substituting (T4"5]) . we obtain 



'■kin 



M 2 Pl 



dtd 3 x Ul + 30 



A - 1 S 8 



+ -(d\d i Bd 2 B + S&BdidjB) 



-{(PhijdkB-ZhijQ&diB 



1 



(^9^)^+0(6 



(5* 



where B is given by (|45|) . and spatial indices are raised and lowered by <5 U and Sij. Note 

that, when written in terms of ( and h^, 

derivatives. 



each term in Ik in includes exactly two time 
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In order to calculate the action up to the (n + 2)-th order (n — 1, 2, • • •), we need 
to solve the momentum constraint up to n-th order. By expanding B and rii as 

B = B 1 + B 2 + --- } n % = nf ] + nf + • • • , (59) 

where B n and n\ n are 0(e n ), and solving the momentum constraint perturbatively, we 
see that B n (and nf) is a sum of various terms with negative powers of (A — 1) up to 
(A — l)~ n (and up to (A — l) - *™ -1 ), respectively) and each term includes just one time 
derivative. This means that Ikin expanded up to 0(e n+2 ) includes various terms with 
negative powers of (A — 1) up to (A — l)~( n+1 ) [jjj and each term includes exactly two time 
derivatives. On the other hand, terms in I 2 =3 i ... i o do not include time derivatives at all 
and are totally independent of A. 

Therefore, while all coefficients of potential terms for ( and h{j remain finite, many 
coefficients of their kinetic terms diverge in the limit A — > 1 + 0. The divergence is 
worse for terms of higher order in the perturbative expansion. This means that the 
naive perturbative expansion breaks down in this limit. Here, let us stress again that 
the theory is still non-perturbatively described by a finite number of parameters, Mpi, 
A, A and Cj (i — 1, 2, • • • , 7) if renormalizable. 

3.4- Non-perturbative continuity at A = 1 + 

Since the naive metric perturbation breaks down in the limit A — > 1 + , nonlinear 
analysis is required. In the following, for simplicity we consider spherically symmetric, 
static, vacuum configurations and show the non-perturbative continuity of the limit. 
In this discussion we consider macroscopic objects and, thus, neglect higher spatial 
derivative terms I z=3 and l z=2 - Anyway, I z=3 and I z=2 have well-behaved perturbative 
expansion and, thus, would not spoil the continuity even if they were included. We set 
the cosmological constant to zero, I z=0 = 0, just for simplicity. 

The lapse is required to be independent of spatial coordinates by the projectability 
condition. Hence, by a space-independent time reparametrization, we can set the lapse 
to unity. Then, by fixing the gauge freedom associated with the spatial diffeomorphism, 
a spherically symmetric, static configuration can be expressed as 

N = 1, Nidx 1 = (3{x)dx, g ij dx l dx j = dx 2 + r(x) 2 dtt 2 2 , (60) 

where dQ 2 is the line element of the unit sphere. The momentum constraint and the 
£:c-component of the dynamical equation are written as 



Br" 

— + (A-1) 



B" B'r' Br" B{r' 






0. 



BB"r 2 + 1B 2 rr" + 4/3/3'rr' + -{B'fr 2 



0, (61) 



1 - (r') 2 + 2BB'rr' + 2B 2 rr" + B 2 {r'Y 

+(A-1) 

where a prime denotes derivative w.r.t. x. The ^-component of the dynamical equation 
follows from the above two equations unless r' = 0, and it is easy to show that r' — 

|| Terms proportional to (A — 1)~(™+ 2 ) cancel after integration by parts. 
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is incompatible with the above two equations for A > 1. We shall not impose the 
global Hamiltonian constraint since we are currently interested in physics in a finite 
region: either staticity or spherical symmetry is not a globally valid assumption and 
thus the equation integrated over a whole space (including e.g. regions far outside the 
cosmological horizon) with these assumptions at face value is not valid. For — 0, the 
second equation leads to r' = ±1 and thus allows only a trivial solution. For this reason, 
hereafter we assume that (3 ^ at least in a neighborhood of a point of interest. 

It is easy to show the continuity of the A — > 1 + limit explicitly. By introducing 
a new variable R(x) by 

R = ^ w )/(») r) ( 62 ) 

we can rewrite equations ( 16"TT) as 

W* (fl ' )2 l -0, (63) 

0. (64) 




A/3 R 

\ /3( A - 1 )A + [(2A-l)/3 2 -l](i? / ) 2 



(65) 



4\R> \(3 ) RR> (3A - l)/3 2 + (A - 1) 

The second equation can be solved w.r.t. /3'//3 and there are two branches: 

^_ 1 ± yi + AAB 
~f3 ~ 2A ' 

(A-l)fl A /3( A - 1 )/ A + [(2A-l)/3 2 -l](iT) 2 

AXR' ' ~ RR' (3A - l)/3 2 + (A - 1) 

The two equations fl63l) and (|65|) provide expressions of highest-order derivatives of R 
and (3, i.e. R" and j3 f , as functions of (R, R', (3). For the '— ' branch, i.e. if we choose 
the ' — ' sign in (165]) . the limit A — > 1 + of the expressions of R and (3 is well-defined as: 



lim R" = 0, lim C= lim - — - ^2 -• (66) 

a^i+o a^i+o /? a^i+o 2i? J R'/3 2 v ; 



These coincide with the equations obtained by simply setting A = 1 in (|63j) and ([61 
Thus, for the '— ' branch, the limit A — > 1 + is continuous. 

For comparison, let us consider general relativity with the metric ansatz 

ds 2 = -dt 2 + [dx + f3(x)dt} 2 + r(x) 2 dQ 2 2 . (67) 

Non-vanishing components of the vacuum Einstein equation G^ = are 

(1- B 2 )(r') 2 - 1 
/V' = 0, W= V W \ (68) 

Remember that we have assumed (3 ^ in a neighborhood of a point of interest. Thus, 
the limit A — > 1 + of the ' — ' branch shown in (16"6"|) agrees with the Einstein equation 
(|68|) . We have thus proved that, for the '— ' branch, the limit A — > 1 + is continuous 
and recovers general relativity for the metric ansatz (1671) . 
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3.5. Schwarzschild solution and Newtonian limit 

The continuity shown above, combined with Birkhoff' s theorem in general relativity, 
implies that the spherically symmetric, static, vacuum solution in the '— ' branch 
approaches a 3 + 1 decomposition of the Schwarzschild spacetime in the A — > 1 + 
limit. This argument neglects higher order spatial curvature terms, I z= 3 and I z =2, but 
this is a fairly good approximation for macroscopic objects. 

If we include I z= % and I z =i then in the A — > 1 + limit we have 

r' = n, J-(^ 2 ) = (r?-l)-E^, (69) 

where r\ is a constant and ct z ir\) (z = 2,3) are constants depending on r\ and the 
parameters in I z=3 and I z =2- Since the spatial metric is flat for r\ = 1, we have 

a 2 (l)=0, = 2,3). (70) 

Integrating (|69|) . we obtain 

r ^ =(r . ? -i )+ ^ + |:gM-L, (7i) 

where /i is an integration constant [23]. For a macroscopic object and thus for large r, 
only the first two terms are important and, as expected, a 3 + 1 decomposition of the 
Schwarzschild spacetime with mass /i is recovered flf 

r' = n, r\{3 2 ~ (r? - 1) + — . (72) 

r 

The 3 + 1 decomposition is characterized by the constant r±. It is noteworthy that for 

n = 1, the solution is not just approximately but exactly the Schwarzschild spacetime 

in the Painleve-Gullstrand coordinate system. This is because the spatial metric is flat 

for ri = 1 and thus higher spatial curvature terms do not contribute to the equations of 

motion (see flTDl ). 

In general relativity the Newtonian limit is usually taken after going to a gauge 
in which the space-dependent part of the lapse is the Newtonian potential. How can 
we express the Newtonian potential in Hofava-Lifshitz gravity with the projectability 
condition? Actually, all information about the Newtonian potential can be included 
in the shift and the spatial metric. See the Schwarzschild solution ( 172|) as an example. 
Even in general relativity, we can choose a gauge in which the lapse is space-independent 
at least locally, and in this gauge the Newtonian potential is encoded in the shift and 
the spatial metric. 

In Hofava-Lifshitz gravity, the same spacetime metric (in the sense of general 
relativity) with different foliations are physically different. Nonetheless, they are 
experimentally and observationally indistinguishable from each other at low energies 
for the following reason. As we all know, Lorentz invariance is a good symmetry of the 

% The Kerr spacetime in a coordinate system with a unit lapse (see e.g. [H]) is also a good approximate 
solution for a macroscopic rotating object. 
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matter sector at least at low energy. It is for this reason that we need a mechanism 
or symmetry to suppress Lorentz violating operators at low energies, as already stated 
at the end of subsubsection 12.1.31 Therefore, although such a mechanism or symmetry 
has not yet been developed and should be explored in detail in the future, we must at 
the very least admit the necessity of recovery of Lorentz symmetry in the matter sector 
at low energy. With this minimal (but challenging) requirement, it is not possible to 
construct low energy observables which can distinguish different foliations of the same 
spacetime through motion of matter. 

In summary, in Hofava-Lifshitz gravity with the projectability condition, the 
Newtonian potential is encoded in the shift and the spatial metric, but matter at low 
energy behaves as if the Newtonian potential were expressed as the space-dependent 
part of the lapse in the "usual" way. Therefore, the projectability condition is not an 
obstacle to expressing the Newtonian potential and taking the Newtonian limit. 

4. Cosmological implications 

There are a number of interesting cosmological implications of Hofava-Lifshitz gravity. 
In this section we shall review some of them: dark matter as an integration "constant" 
(subsection 14. ip . bouncing and cyclic universes (subsection I4.2[) and generation of scale- 
invariant cosmological perturbation from z = 3 scaling (subsection 14. 3p . 

4-1. Dark matter as an integration "constant" 

4-1.1. Structure of GR and FRW universe General relativity has the four-dimensional 
diffeomorphism invariance as its fundamental symmetry. As a result, there are four 
local constraints: one Hamiltonian constraint and three momentum constraints at each 
spatial point and at each time. The constraints are preserved by dynamical equations. 
Thus, we can solve dynamical equations without worrying about constraints, provided 
that constraints are satisfied at initial time. 

Now, let us consider the flat FRW spacetime 

ds 2 = -dt 2 + a 2 (t)dx 2 , (73) 

as a simple example. It is supposed that this metric approximates overall behavior of 
our patch of the universe inside the Hubble horizon. The Hamiltonian constraint leads 
to the Friedmann equation 

3H 2 = 87rGp, (74) 

where G is Newton's constant, H = a/a is the Hubble expansion rate and p is the total 
energy density of matter contents of the universe. Equations of motion of matter lead 
to the conservation equation 

p + 3H(p + P) = 0, (75) 
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where P is the total pressure of the matter contents. The momentum constraint is 
trivial because of the symmetry of the FRW spacetime. The dynamical equation 

- (2H + 3H 2 ) = 8nGP (76) 

follows from the Friedmann equation and the conservation equation, and thus we do not 
consider it as an independent equation. 

4-1-2. Structure of HL gravity and FRW universe The fundamental symmetry of 
Hofava-Lifshitz gravity is the invariance under the foliation preserving diffeomorphism 
(jlip . which is 3-dimensional spatial diffeomorphism plus space- independent time 
reparametrization. Consequently, contrary to general relativity, the theory has 3 local 
constraints and 1 global constraint: 3 momentum constraints at each spatial point at 
each time and 1 Hamiltonian constraint integrated over a whole space at each time. 
Of course, the constraints are preserved by dynamical equations. Thus, we can solve 
dynamical equations without worrying about constraints, provided that constraints are 
satisfied at initial time. 

Now let us consider the flat FRW spacetime (1731) . or 

N = l, N l = 0, 9ij = a(t) 2 6ij, (77) 

in Hofava-Lifshitz gravity. Again, the FRW spacetime is supposed to approximate 

overall behavior of our patch of the universe inside the Hubble horizon. This means 

that the global Hamiltonian constraint, which is an integral over a whole space including 

regions far outside the Hubble horizon, does not apply to our system within the horizon. 

Thus, the lack of local Hamiltonian constraint implies that there is no Friedmann 

equation and that we should consider the dynamical equation 

3A — 1 
-— — - (2H + 3H 2 ) = 8tiGP (78) 

as an independent equation. Here, note that higher spatial curvature terms do not 
contribute to the equation because of the spatial flatness. Equations of motion for 
matter leads to the conservation equation (I7"5l) at least at low energy, provided that 
the local Lorentz invariance is restored in the matter sector at low energy as required 
by many experimental and observational data (see discussion in the second-to-the-last 
paragraph of subsection I3.5p . At high energy, however, the matter sector does not have 
to satisfy the conservation equation and thus the equation of motion for matter generally 
leads to 

p + 3H{p + P) = -Q, (79) 

where Q represents the amount of energy non-conservation. Note that Q — > at low 
energy. The two equations ( 1751) and ( 1791) are sufficient to describe the evolution of our 
system. Indeed, it is easy to obtain the first integral of the dynamical equation: 

C{t) 



3 ( 3A -V = 8,G 



P + 



a 3 



(80) 
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where 

C(t) = C + f Q(t')a 3 (t')dt', (81) 

J to 

and C = C(t ) is an integration constant. Since Q — > at low energy, 

C(t) — >• const, at low energy. (82) 

The first integral ( 1HU1) looks like Friedmann equation but, intriguingly, the extra term 
(oc C(t)/a 3 ) behaves like dark matter at low energy. This term is not real matter 
but gravitationally behaves like pressure- less dust. Thus, in Hofava-Lifshitz gravity, 
something like dark matter emerges as an integration "constant" at least in flat FRW 
background at low energy. Note that (jHTl) describes how the "dark matter" is generated 
in the early universe: even with C$ = and to — — oo, we have non-vanishing C(t) at 
late time. 

4-1.3. General case in the IR We now show that the dark matter as an integration 
"constant" emerges at low energy in more general situation. 

Low energy behavior of the theory is described by the IR action (123]) . This looks 
like the Einstein Hilbert action with the ADM decomposition if A = 1. Hence, we set 
A = 1 in the discussion below, hoping that in the near future we can show that A = 1 is 
a stable IR fixed point of the RG flow. 

The Hamiltonian constraint is then of the form 

J Sx^g [G$ + A<7$ - SnGT^} n»n v = 0, (83) 

where gjfj is the four-dimensional spacetime metric defined in ( Tl3l) . G^J is the Einstein 
tensor of gffl, n M is the unit vector normal to the constant time hypersurface defined in 
(133J1 . and T^ v is the energy momentum tensor of matter. This is an equation integrated 
over a whole space including regions far outside the cosmological horizon, and thus does 
not restrict physics inside our patch of the universe. On the other hand, the momentum 
constraint 

0, (84) 

(85) 

are local equations. 

Interestingly, it is possible to give a general solution to these local equations. For 
this purpose let us define T^ rk by 

T^ k = (SttG)- 1 [G{S + Ag$ - 8ttGT, u ] . (86) 

The momentum constraint (184"|) and the dynamical equations fl85l) are rewritten as 

T dark n » = 0> T dark = q^ ( g? ) 



19 





G { $ + Ag$ - MGT Hl _ 


n^ 


and dynamical equations 


( 


lf)+kgf-&<KGT ij = 


-o, 



meaning that the time-space components and the space-space components of T^ rk 
vanish. Only the time-time component remains and thus T^ rk should be proportional 

to n^riy: 

TpJ = Pdark r >' l _ i n u , (88) 

where the scalar pd ar k can in general depend on both time and spatial coordinates. This 
is exactly of the form of pressure- less dust and thus behaves like dark matter. It is easy 
to show that the vector n 1 * defined in ( 133]) follows the geodesic equation 

n u V v n ,M = 0. (89) 

Also, by taking the divergence of the definition ( 186]) of T^ rk , we can show that Pdark 
satisfies the conservation equation 

n^d^Pdark + K Pdark = 0, (90) 

provided that the real matter sector recovers the local Lorentz invariance in the IR 
and thus satisfies the energy conservation n^ v T ilv = at low energy. In more general 
cases the right hand side of ( 1901 obtains non-vanishing contributions from higher spatial 
curvature terms, deviation of A from unity, and energy non-conservation of matter. 

As a consistency check, let us apply the conservation equation ( 1901) to the flat 
FRW spacetime (1771) . In this case, (190]) is reduced to d t pdark + 3Hp dark = and thus 
Pdark oc a~ 3 . This reproduces the scale factor dependence of the last term in ( !80l) with 
C(t) = const. 

In summary, we have shown that gravity equations of motion in Hofava-Lifshitz 
gravity at low energy with A = 1 is written as 

G$) + A($ = 87rG % v + PdarkU^} . (91) 

This modified Einstein equation includes a built-in component which behaves like dark 
matter, as an inevitable consequence of the projectability condition. The "dark matter 
velocity vector" n M follows the geodesic equation (18~9~1) and the "dark matter energy 
density" pdark satisfies the conservation equation (190]) . In the Newtonian limit the 
modified Einstein equation (19T|) reduces to the Poisson equation with the built-in "dark 



matter" included. Note that, as already discussed in subsection 13.5] the Newtonian 
potential is encoded not in the lapse but in the shift and the spatial metric. 

4-2. Bouncing and cyclic universes 

Higher curvature terms in the action are expected to play important roles in the early 
universe. In this section we consider the FRW universe with spatial curvature 

dr 2 
l-Kr 2 

and see that higher curvature terms drastically change the evolution of the early 
universe. In particular, bouncing universes and cyclic universes are allowed as regular 
solutions in Hofava-Lifshitz gravity. 
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N = 1, N { = 0, g ij dx i dx j = a(tf 



+ r 2 dn 2 9 



(92) 



4-2.1. Modified Friedmann equation with higher curvature terms As already stated in 
subsubsection I4.1.2[ the FRW spacetime is just an approximation to describe overall 
behavior of our patch of the universe inside the Hubble horizon, and thus the global 
Hamiltonian constraint does not restrict the dynamics of our approximate FRW universe 
inside the horizon. (This is evident in the presence of superhorizon fluctuations.) 
Instead, we should consider the dynamical equation as an independent equation, 

- 3 -± y ± (2H + 3H i ) = S ,GP-^-^ + !i-A, (93) 

2 a b a 4 a z 

where 

a 3 = 1927rG(c 3 + 3c4 + 9c 5 ), a 2 = 32tiG{c 6 + 3c 7 ). (94) 

By using the definition (179|) of energy non-conservation Q, we can easily obtain the 
first integral of the dynamical equation, 

] -H 2 = 8ttG 



3(3A-1)„ 2 



C{t) 



"» K * 3a * K2 3K + A, (95) 



a 6 a 4 a 2 



2 

where C(t) is defined in (181]). This is a straightforward generalization of ( 180]) and 
includes contributions from the spatial curvature K. As in ( IHUj) . the term proportional 
to C(t)/a 3 behaves like dark matter at low energy as C(t) — > const. In the early 
universe, i.e. for small a, the curvature cubic term (oc K 3 /a 6 ) plays important roles. 

In order to see qualitative behavior of the system, let us rewrite the first integral 
( 195]) in the form of the energy conservation equation for a non-relativistic particle moving 
in a 1-dimensional potential as 

l -d 2 + ^jV(a) = 0, (96) 



where 



„ f v a 3 ^ 3 , a 2 K 2 K A 2 4vrG 

via) = — = — a 

v ; 6a 4 2a 2 2 6 3 



2 C(t) 
a 



(97) 



The shape of the potential V(a) completely determines the behavior of the system. 

4-2.2. Simple examples Let us now consider some simple examples. For simplicity we 
set «3 = 1, a% — 0, K — 1, p — 0, C — const. We still have freedom to choose values 
of A and C. We show four examples of the 1-dimensional potential ( 197]) : a bouncing 
universe (Figured]), a cyclic universe (Figure [2}, an unstable static universe (Figure [3]) 
and a stable static universe (Figure H]). See [25] for more examples with p = and 
C = const. 

4-3. Scale-invariant cosmological perturbations from z = 3 scaling 

One of the essential ingredients of Hof ava-Lifshitz gravity is the anisotropic scaling with 
the dynamical critical exponent z > 3. Indeed, it is this property that makes the theory 
power-counting renormalizable and attractive as a candidate for the theory of quantum 
gravity. There are interesting cosmological implications of the anisotropic scaling. In 
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Figure 1. V(a) for 

A = 0.4 and 8ttGC = 
10. If the universe is 
initially contracting then 
it bounces and expands. 



Figure 2. V(a) for A = 
and 8ttGC = 10. This 
allows a series of expansion 
and contraction, i.e. a 
cyclic universe. 





Figure 3. V(a) for 

A = 15/64 and SttGC = 
17/4. This allows a static 
universe at a — 2 but it 
is unstable. A bouncing 
universe is also allowed. 



Figure 4. V(a) for A = 
and C = 4. This allows 
a stable static universe at 

= 1. 
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this section we show that the anisotropic scaling with the minimal z, i.e. z = 3, leads 
to a new mechanism for the generation of scale-invariant cosmological perturbations. 
Intriguingly, this mechanism works even without inflation. 

4-3.1. Usual story with z = 1 Before explaining the new mechanism, let us remind 
ourselves of the usual story with z = 1. 

Cosmological perturbations are analyzed by perturbative expansion around a FRW 
background. In the linearized level, perturbations are Fourier expanded and the 
evolution of each mode is characterized by the frequency u defined by the dispersion 
relation 

u 2 = c 2 ^ (98) 

where c s is the sound speed, k c is the comoving wave number and a is the scale factor 
of the universe. For simplicity we assume that the time dependence of c s , if any, is 
slow compared with the cosmological time scale H~ x , where H = a/ a is the Hubble 
expansion rate. (For example, c s is identically 1 for a canonical scalar field with any 
potential.) 

If a mode of interest satisfies u 2 ^> H 2 then the evolution of the mode is not affected 
by the expansion of the universe and the mode just oscillates. When u 2 <C H 2 , on the 
other hand, the expansion of the universe is so rapid that the Hubble friction freezes the 
mode and the mode stays almost constant. Generation of cosmological perturbations 
from quantum fluctuations is nothing but the oscillation followed by the freeze-out. 
Therefore, the condition for generation of cosmological perturbations is 

Jt b) > °" (99) 

With the z — 1 dispersion relation ( |98|) . this condition is equivalent to a > for 
expanding universe (d > 0). Therefore, if z = 1 then generation of cosmological 
perturbations from quantum fluctuations requires accelerated expansion of the universe, 
i.e. inflation. For example, for power law expansion a oc t p , p > 1 is required. 

Observational data of the cosmic microwave background strongly indicates that 
the primordial cosmological perturbations have an almost scale-invariant spectrum. It 
is easy to see that the scale-invariance also requires inflation. From the scaling fl2]) with 
s — 1, the amplitude of quantum fluctuations of the scalar field should be proportional 
to the energy scale of the system. In cosmology the energy scale is set by the Hubble 
expansion rate H. Thus, we expect that 

5<f) oc H. (100) 

Since cosmological perturbations with different scales are generated at different times, 
the scale-invariance is nothing but the constancy of the right hand side of fllOOp . Noting 
that H = a/ a, this implies the exponential expansion of the universe a oc exp(Ht), 
namely inflation. 
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We have seen that, for z = 1, both the generation of cosmological perturbations and 
the scale-invariance of generated perturbations require the existence of an inflationary 
epoch in the early universe. 

4-3.2. The story in the UV with z = 3 The condition ( 1991) for generation of cosmological 
perturbations is valid irrespective of the dispersion relation. In Hof ava-Lifshitz gravity, 
to realize the anisotropic scaling ([5]), the dispersion relation for a physical degree of 
freedom in the UV should be 

"'= m2 *(M' (ioi) 

where M is some energy scale. By substituting this to the condition ( 199]) we obtain 
d 2 (a z )/dt 2 > for expanding universe (a > 0). Since z > 3 in Hofava-Lifshitz gravity at 
high energy, generation of cosmological perturbations from quantum fluctuations does 
not require accelerated expansion of the universe, i.e. inflation. For example, power law 
expansion a oc t p with p > 1/z satisfies the condition. 

In this way, the anisotropic scaling provides a solution to the horizon problem. 
Essential reason for this is that perturbations freeze-out not at the Hubble horizon but 
at the sound horizon, defined by u ~ H . The physical radius of sound horizon is thus 
~ (M z ~ l H)~ l l z . In the UV epoch (H ^> M), the sound horizon is far outside the 
Hubble horizon and can therefore accommodate scales much longer than the Hubble 
horizon size. In order to stretch microscopic scales to cosmological scales, we just need 
to have a long enough expansion history (satisfying the condition d 2 (a z )/dt 2 > 0) in 
the UV epoch. Note that M is not a cutoff scale of a low energy effective theory but 
is just the scale at which the theory starts exhibiting the anisotropic scaling, provided 
that Hofava-Lifshitz gravity is UV complete. 

For general z, the formula (JTj) implies that the amplitude of quantum fluctuations 
of <p should be 

3-z 

5(j) ^Mx^y\ (102) 

where M is defined through the dispersion relation f llOip . This is of course consistent 
with the well-known result (HOOp for z = 1 and the result in ghost inflation 5<p ~ 
(M 3 H ) 1 / 4 [261 I2Z] for z — 2. On the other hand, in Hofava-Lifshitz gravity with the 
minimal value of z ,i.e. z = 3, f l 1 2 f) is reduced to 5<p ~ M, implying that the amplitude 
of quantum fluctuations does not depend on the Hubble expansion rate. This means 
that the spectrum of cosmological perturbations in Hofava-Lifshitz gravity with z = 3 
is automatically scale-invariant even without inflation. 

4-3.3. A simple model We have shown that the anisotropic scaling with z = 3 naturally 
leads to a new mechanism for generation of scale-invariant cosmological perturbations. 
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As a simple implementation of the mechanism, let us consider a free scalar field described 
by the action 



- / dtd A xN^fg 



-L^-wd^y + tpotp 



(103) 



where 



° = §i-^ + c >-< A = 9 lJD ^ (104) 

This is a covariantized version of (JHJ). 

In the UV, the first term in O is dominant and the scalar field action exhibits 
the z = 3 scaling. In this regime it is easy to find the mode function in a flat FRW 
background as [5] 

x 2-V^-^Mexp l-i^ % , (105) 



(*0Vtf=(|) • ( 106 ) 



where a is the scale factor, t is the proper time, k c is the comoving wave number 
and k c = \k c \. Note that this is not just WKB approximation but actually exact and 
applicable to both subhorizon and superhorizon scales in any background a(t), provided 
that the first term in O is dominant. The mode function approaches a constant value in 
the a — > oo limit if and only if the integral J°° dt/a 3 converges. For power-law expansion 
a oc t p , this condition is satisfied if p > 1/3, and agrees with the condition for the freeze- 
out after oscillation discussed after fllOll) . Provided that the integral converges, the 
power-spectrum is calculated as 

* 27T 2 

This is manifestly scale- invariant in accord with the general argument after ( J102R . In 
this way, scale-invariant cosmological perturbations of the scalar field can be generated 
even without inflation. 

After scales of interest exit the sound horizon, cosmological perturbations of the 
scalar field can be converted to curvature perturbations by either curvaton mechanism or 
modulated decay of heavy particles or/and oscillating fields. For example, it is possible 
to suppose that the scalar field itself plays the role of a curvaton [5] . When the Hubble 
expansion rate becomes as low as m<p, <fi starts rolling and eventually decays to radiation. 
Perturbations of <fi are converted to those of radiation energy density and thus curvature 
perturbations. 

In the IR, the first two terms in O can be neglected and the usual z = 1 scaling is 
recovered. In this epoch, unless the universe is in an inflationary phase, physical scales 
re-enter the horizon as usual. 

5. Summary and discussions 

We have reviewed basic construction and cosmological implications of a power-counting 
renormalizable theory of gravitation recently proposed by Hofava. While there are 
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many fundamental issues to be addressed in the future, it is interesting to investigate 
cosmological implications. 

Since the high energy behavior of Hofava-Lifshitz gravity is very different from 
general relativity, there is a possibility that the theory does not exactly recover general 
relativity at low energy. As reviewed in subsection 14.11 this is indeed the case and the 
theory can instead mimic general relativity plus dark matter. The constraint algebra in 
this theory is smaller than general relativity since the time slicing is synchronized with 
the "dark matter rest frame" in the theory level. In subsection 14.21 we have shown that 
higher spatial curvature terms in the action drastically change the evolution of the early 
universe. We have derived modified Friedmann equation with higher spatial curvature 
terms and have shown some simple examples, including bouncing and cyclic universes. 
The anisotropic scaling at high energy is one of essential ingredients of the theory since 
the power-counting renormalizability stems from it. In subsection 14.31 we have reviewed 
a new mechanism for generation of cosmological perturbations based on the anisotropic 
scaling. This mechanism can solve the horizon problem and generate scale-invariant 
cosmological perturbations even without inflation. 

In Sec. [3] we have commented on some issues related to the scalar graviton and the 
A — > 1 + limit, where A is a parameter in the kinetic action. We have explicitly seen 
that the naive metric perturbation breaks down for the scalar graviton in the A — > 1 + 
limit. However, this does not necessarily imply the loss of predictability. Actually, for 
spherically-symmetric, static, vacuum configurations we have proved that the limit is 
non-perturbatively continuous and safely recovers general relativity. 

Now let us compile a list of some important open questions. 

• Renormalizability must be shown beyond power-counting argument. (See [28] 
for discussion about renormalizability of the theory with the detailed balance 
condition.) 

• The RG flow of the theory must be analyzed. In particular, it is very important to 
see whether A = 1 is an IR fixed point or not. If it is the case then we would like 
to know whether the RG flow can satisfy the condition (|56|) or not. 

• We have to develop mechanisms or symmetries to suppress Lorentz violating 
operators in the matter sector at low energies. Perhaps, embedding into a larger 
theory is needed. One such possibility is related to supersymmetry [16J. 



Is there an analogue of Vainshtein effect [H]? In subsection 13.41 non-perturbative 
continuity of the A — > 1 + limit was shown only for spherically-symmetric, static, 
vacuum configurations. We need to consider more general situations in order to see 
how general the non-perturbative continuity is. 

In [10], based on exact results in some simple cases, it was conjectured that there 
is no caustic for constant time hypersurfaces. We need to provide evidences for 
this conjecture in more general situations if a proof is difficult. Perhaps, numerical 
simulations similar to those in [29] are necessary. 
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• In |23j it was proved that a spherically-symmetric solution should include a time- 
dependent region near the center. On the other hand, as shown in Sec. [3] of 
the present article, the vacuum region far from the center recovers the standard 
Schwarzschild geometry. Since the size of the dynamical region is expected to 
be of the fundamental scale, the dynamical nature of the central region is not 
really relevant for macroscopic objects such as astrophysical stars. Microscopically, 
however, this could be rather significant . We would like to know, e.g. the typical 
size of the dynamical region and the motion of its boundary. 

• As already stressed in [TO] , we need to know if microscopic lumps of "dark matter 
as an integration constant" can play the role of particles in usual dark matter 
models from macroscopic viewpoint. Interactions among them such as collisions 
and bounces need to be understood. At astrophysical scales, we need to see if 
collective behavior of a group of large number of microscopic lumps can more or 
less mimic behavior of a cluster of particles with velocity dispersion and vorticity. 
Clearly, detailed investigation is necessary to understand rich dynamics of "dark 
matter" from microscopic to macroscopic scales. 

• As shown in IJ8T1) . "dark matter as an integration constant" is generated in 
the early universe even if it vanishes initially. This formula can be applied to 
superhorizon perturbations. Given a concrete model of the matter sector, therefore, 
it is straightforward to estimate the typical amplitude and spectrum of the "dark 
matter" . If a single physical degree of freedom is responsible for both the source 
term Q in ( 1HT1) and generation of curvature perturbations then it should be possible 
to realize adiabatic initial conditions for the late time evolution of perturbations. 
(See [30J for classical late time evolution.) It is worthwhile investigating this 
possibility in details. 

• The mechanism reviewed in subsection 14.31 generates scale-invariant cosmological 
perturbations without a need for inflation. It would be interesting to see whether 
renormalization effects such as anomalous dimension can break the exact scale- 
invariance and explain the observed spectral tilt. 

• In the early universe, it is expected that A should deviate from 1 under the RG flow 
and that the scalar graviton can be treated perturbatively. Since the scalar graviton 
has the z = 3 anisotropic scaling in the UV, it should also obtain the scale-invariant 
cosmological perturbations |3~T] . Provided that A = 1 is a stable IR fixed point of 
the RG flow, as the universe expands and the Hubble expansion rate decreases, 
A approaches 1 and the perturbative treatment of the scalar graviton becomes 
invalid □. However, the result in subsection 13.41 suggests that the A — > 1 + limit 
may be non-perturbatively continuous. A natural question is then "how to convert 
the scale-invariant cosmological perturbations of the scalar graviton to observables 

+ Thus, the conventional cosmological perturbation scheme |32j probably breaks down for the scalar 
graviton in the A — > 1 + limit. Again, this does not necessarily imply loss of predictability but requires 
nonlinear analysis. 
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such as cosmic microwave background anisotropics and matter power spectrum?" 

Since Hofava's original proposal in January 2009, several extensions appeared in the 
literature. Bias, et.al. [21] proposed an extension without the projectability condition by 
including spatial derivatives of the lapse in the action. More recent proposal by Hofava 
and Melby-Thompson [33J respects the projectability condition but the fundamental 
symmetry of the theory is larger than the original one. 

Throughout this article, we have considered the minimal theory, i.e. the original 
theory with the projectability condition but without extension of the symmetry. 
Whether this minimal theory is viable is still an open question and crucially depends 
on non-perturbative nature of the scalar graviton (see subsection 13.41) and properties of 
the RG flow (see the condition (1561) ). 
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